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Abstract. In specific types of partially rectangular billiards we estimate the mass of 
an eigenfunction of energy E in the region outside the rectangular set in the high-energy 
limit. We use the adiabatic ansatz to compare the Dirichlet energy form with a second 
quadratic form for which separation of variables applies. This allows us to use sharp one- 
dimensional control estimates and to derive the bound assuming that E is not resonating 
with the Dirichlet spectrum of the rectangular part. 



1. Introduction 

We study concentration and non-concentration of eigenfunctions of the Laplace opera- 
tor in stadium-like billiards. As predicted by the quantum/classical correspondence, such 
concentration is deeply linked with the classical underlying dynamics. In particular, the 
celebrated quantum ergodicity theorem roughly states that when the corresponding clas- 
sical dynamics is ergodic then almost every sequence of eigenfunctions equidistributes in 
the high energy limit (see [151 S HE] and [HI [17] in the billiard setting for a more precise 
statement). In strongly chaotic systems such as negatively curved manifolds, it is expected 
that every sequence of eigenfunctions equidistributes. This statement is the Quantum 
Unique Ergodicity conjecture (Q.U.E.) and remains open in most cases despite several re- 
cent striking results (see for instance [H [121 E l2])- On the other extreme, the Bunimovich 
stadium, although ergodic, is expected to violate Q.U.E.. Indeed, it is expected that there 
exist bouncing ball modes i.e. exceptional sequences of eigenfunctions concentrating on the 
cylinder of bouncing ball periodic orbits that sweep out the rectangular region (see [5] for 
instance). The existence of such bouncing ball modes is still open and only recently did 
Hassell prove that the generic Bunimovich stadium billiard indeed fails to be Q.U.E. (see 

nun. 



Our work is closely related to the search for bouncing ball modes but proceeds loosely 
speaking in the other direction. We actually aim at understanding how strong concen- 
tration of eigenfunctions in the rectangular part cannot be. We thus follow [7J in which 
Burq-Zworski proved that even bouncing ball modes couldn't concentrate strictly inside 
the rectangular region. This was made precise by Burq-Hassell-Wunsch in [B] where the 
following estimate was proved : 

IMIl2(wq > -E^IMIl 2 ^) 
l 
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in which ||m||l 2 (vk) (resp. ||w||z,2(f2)) denotes the L 2 norm of the eigenfunction u in the wings 
(resp. in the billiard). 

Our main result for the Bunimovich stadium is the following theorem 

Theorem 1. Let Q be a Bunimovich stadium with rectangular part R := [— Bq, 0] x [0, Lo]. 
We set W = Q \ R and denote by E the Dirichlet spectrum of R, i. e. 

k 2 7T 2 l 2 TX 2 



T2 1 D2 ' 

^0 n o 

For any e > there exists Eq and C such that if u is an eigenfunction of energy E such 
that E > Eq and dist(i?, E) > E~ e then the following estimates holds: 

IM|l 2 (X!) <CE^\\u\\ L 2 {w)) 

This bound improves on the Burq-Hassell-Wunsch bound provided that e < |. It is 
natural that the smaller e is the better the bound is. Indeed, the condition on the distance 
between E and E is comparable to a non-resonance condition and should imply heuris- 
tically that u must have some mass in the wing region. It is quite interesting to have 
a quantitative statement confirming this heuristics. We will actually give a more general 
statement concerning more general billiards (see Theorem [2]). In particular we will consider 
billiards with smoother boundaries (see Def. E]) disregarding the fact that these may not be 
ergodic. Here again we expect the bound to be better when the billiard becomes smoother 
and this statement is made quantitative in Theorem [2J 

The method we propose relies on comparing the Dirichlet energy quadratic form with 
another quadratic form arising from the adiabatic ansatz presented in the numerical study 
of eigenfunctions by Backer-Schubert-Stifter [3]. This adiabatic quadratic form has also 
appeared recently in the works of Hillairet- Judge [11] in the study of the spectrum of the 
Laplacian on triangles. These two quadratic forms are close provided we do not enter too 
deeply into the wing region so that the non-concentration estimate really takes place in 
a neighbourhood of the rectangle that becomes smaller and smaller when the energy goes 
to infinity (see Sections 14.3.31 and 14.6. ip . Since the new quadratic form may be addressed 
using separation of variables, we will show precise one-dimensional control estimates and 
then use them to prove our results. We have separated these one dimensional estimates 
in an appendix since they may be of independent interest. Finally, we remark that the 
method can be applied to quasimodes with some caution (see Remark 15. 2p but there are 
no reasons to think that the bound we obtain is optimal. 
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2. The setting 

Take L a function defined on [—-Bo, B\\ with the following properties : 

- For non-positive x, L(x) = L > 0. 

- On (0,Bi), L is smooth, non-negative and non- increasing. 

- When x goes to Bi, V has a negative limit (either finite or — oo). 

- For small positive x, we have the following asymptotic expansions: 

, f L(x) = L - c L x 7 + o(x 7 ), 

' \ L'(x) = -CL-fX^' 1 + o(x^- 1 ) 

for some positive cl and 7 > 3/2. 
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Figure 1. An example of the billiard fl. 

The billiard Q is then defined by 

Q = {(x,y), I - B < x < B u < y < L(x) } . 

See Fig. [T]for an example of an applicable billiard. For any b < B x , we will denote by 
Q b :=nn{x < b} and by W b := SI n {0 < x < b}. 
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We study eigenf unctions of the positive Dirichlet Laplacian, A, on Q. Namely, we study 
solutions, ue such that 

Amb = - {pi + dy) u E = Eu E , 

uslan = 0, 

where E > 0. 

We may formulate this equation using quadratic forms. We thus introduce q defined on 
H\Q) by 

q{u) = I | Vu\ 2 dxdy. 
Jn 

The Euclidean Laplacian with Dirichlet boundary condition in Q is the unique self-adjoint 
operator associated with q defined on Hq(Q). We denote by q b the restriction of q to if ^f^) 
and by A;, the Dirichlet Laplace operator on We will also denote by V b the set of smooth 
functions with compact support in fl b - 



3. Adiabatic approximation 



Motivated by the well-known eigenvalue problem on a rectangular billiard and compu- 
tational results in j3], we introduce a second family of quadratic forms a b and compare it 
to q b . 

For any b < B\ and any u G V b , Fourier decomposition in y implies that 



(3.1) 
Since 



u[x 



y) = ^Mfc(x) sin 

k 



L{x) 



sin ( kn 



L(x) 



nk 



My 



y 



L(x) 



2 




[L{x)_ 


I 



each Fourier coefficient u k is given by 

u k (x) = 
For such u, we define 



u(x, y) sin 



uk 
Ux) 



V dy. 



2 kW 



-B 



u' k {x)\> + 



L 2 (x) 



\u k (x)\' 



L(x) 



dx, 



L(x] 



\u k (x)\' 



dx. 
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Observe that for each fixed x, Plancherel's formula reads 

EK(x)| 2 ^ = [ LiX) \u(x,y)\ 2 dy, 
fceN ^° 

so that we get iV&(«) = || w||| /2 ^ b ) by integration with respect to x. 

Fixing some < bo < B\, and using that L is uniformly bounded above and below on 
[—Bq, bo] we find a constant C such that for any b < bo and u G L 2 (Qb) : 



(3-2) C" 1 ||«ll^<ElKHi 2 (-Bo ! 6)< C 'll u 



2 

ir"fc|lL 2 (-B ,b) — ^ H"ll^( 
fc=l 



The quadratic form a& appears as the direct sum of the following quadratic forms a^k 
(that can be defined on the whole function space H 1 ^— B , b)) : 

(3.3) a b;k (u):= (\uf + — ^{uf ) -^-dx. 



Bo 



L 2 (x) 



Recall that, on an interval J, the standard H l norm is defined by 



(3-4) \\u\\ H i := [\\u'\\ 2 L2{I) + ||w|||2 (J) 

so that, for any k and b < B\ and any u G B , b) we have 

/ k 2 TX 2 \ / k 2 7T 2 \ 

(3.5) min^L(6), J \\ u \\m < a &,*W < max , IMIffi- 

i 

The norm thus defines on H 1 (—B ,b) a norm that is equivalent to the standard H 1 
norm. 

3.1. Comparing a& and g&. To compare a& and g&, we introduce the following operators 
.D and R defined on T>^ by 

ita = 

Du = d x u + 

Using Plancherel formula for each fixed x and then integrating, we obtain 

ctb(ii) = / \Du\ 2 + \d y u\ 2 dxdy . 
from which the following holds for any u, v G T>b. 
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a b (u,v) - q b (u,v) = (Du, Dv) - (d x u, d x v) 

(3.6) = (d x u, Rv) + (Ru, Dv) , 

(3.7) = (d x u,Rv) + (Ru,d x v) + (Ru,Rv). 

We thus obtain the following lemma. 
Lemma 3.1. Let 5 be the function defined by 

5(b) = sup|L'(x)| + sup|Z/(x)| 2 . 

(0,6] (0,6] 

Then for all u, v G V b 

\a b (u,v) - q b (u, v)\ < 5(b) ■ q b 2 (u) ■ q b 2 (v). 

Remark 3.1. The function 5 is continuous on (0,Bi) and 5(b) = 0(6 7_1 ) when b goes to 
0. 

i 

Proof. In (13 .7p . we use Cauchy-Schwarz inequality, max(||Du||, ||9j/u||) < a b (u), and the 
fact that is uniformly bounded by 1 on Q. □ 

The following corollary is then straightforward. 

Corollary 3.2. For any < b < B\ and any u G if 1 (fi), the linear functional A defined 
by A(v) := a b (u,v) — q b (u,v) belongs to H' 1 ^^. Moreover 

||A||#-i(n 6 ) < 5(b)\\u\\ H i(n b y 

4. Non- Concentration 

4.1. Preliminary reduction. Let u be an eigenfunction of q with eigenvalue E. And 
define the associated linear functional A using corollary 13.21 

Integration by parts shows that for any v G ifg(fi&) we have 

q b (u,v) = E ■ (u,v) L2{n) , 

so that we have 

(4.1) a b (u,v) - E ■ N b (u,v) = A(v). 

We now deal with this equation using the adiabatic decomposition. We thus define A^ 
as the distribution over T> b such that, for any v G D bl 

(4.2) A k (v) := A ( v(x)sin I kn 



L(x) 
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Remark 4.1. From now on, u will always denote the eigenfunction that we are dealing 
with. We will denote by u k the functions entering in the adiabatic decomposition of u, by 
A the linear functional associated with u and by A k the one-dimensional linear functionals 
that are associated with A. 

A straightforward computation yields, that for any v EV b we have 

f b L(x) 
a h)k (u k ,v) - E ■ j u k (x)v(x)——dx = A k (v), 

where a b)k is the quadratic form defined in fl 3 . 3 [) . 

An integration by parts then shows that, in the distributional sense in (—Bo, b), we have 

(4.3) -__ (K)+ ^_- £ j Iit = Ai , 
where the linear functional A k is defined by 

(4.4) A k (v) :=A k (l^y 

Remark 4.2. Since L is not smooth, this definition of A^ doesn't make sense as a distri- 
bution. However, in the next section, we will prove that A k actually is in if -1 and, since 
multiplication by j is a bounded operator from i/ _1 (— B ,b) into itself, we thus get that 
A k is a perfectly legitimate element of H" 1 . Moreover, for any fe there exists C(b ) such 
that for any b < bo, and v G T> b , we have 

2 

||-^||tf-H-Bo,&) < C(bo)\\v\\ H -i { _ B(hb) . 
We denote by P k the operator that is defined by 

n / N 1 d ,r /N ( k ^ 2 t\ 

fl b („) = - I -(L,0 + ^_-^ Ul 

and we try to analyze the way a solution to equation (14. 3D on (—Bo, b) may be controlled 
by its behaviour on (0, b). 

The strategy will depend upon whether k is large or not, but first we have to get a bound 
on A k in some reasonable functional space of distributions. 

4.2. Bounding A k . In this section, we prove that each A*, is actually in H~ 1 (—B , b) and 
provide a bound for its H" 1 norm. 

We first note that, using (ET4"1) . for any F G H~ l (-B ,b) : 

\F((j>)\ IW)I 

(4-5) WFWH-H-Bofi) ■= sup < sup \-£r > 



4>ev b \\9WW- <t>&v b W<P Hl 2 
Using (13. 6 p in the definition of A k -see (14. 2ft - we obtain 



s 
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A k (v) = (d x u, R ( v(x) sin ( kir Jj ) ) ) + (Ru, D ( v(x) sin ( kir 



L x) 



L{x) 



Denote by A k (v) the first term on the right and B k {v) the second term. By inspection, 
we have 



where we have set 
(4.6) F k (x) 



A k (v) '■— [ v{x)^^-F k (x)dx 
B k (v) :=~J v'(x)L'(x)G k (x)dx, 



(4.7) 



G k (x) 



f {x > ( y \ 

j l w ■ yd x u(x, y) ■ cos Ik— r J dy 
L{x)l lw - ydA ^ v) -^{^W)) dy 



L(x 
2 r L ( x ) 



Since u G if 1 (f2), F k and G k are L 2 (0, b) and we can estimate the H 1 norm of A k using 
them. 

Lemma 4.1. For any bo < B\, and given A k and F k , G k defined as above, there exists 
C = C(Q b0 ) 

(4.8) ||A fc || H -i < C(W\\F k \\ L2m + b^WGkW 

Proof. We estimate A k (v), using first an integration by parts 



a / \ kir 
A k v) := — - 



v \x) 



B 



HO 



F fc (0 d£ dx. 



Using Cauchy-Schwarz inequality, and the fact that L'(£) = 0(£ 7 x ) we have 



L'(0 



F k (Od£ 



< Cx\ a ||F fc ||i2(o l6 ). 



/o m 

Inserting into A k (v) and using Cauchy-Schwarz inequality again we get 

\A k (v)\ < C ■ (A;6 7 )||F fc || L2(0i6) • \\v'\\ L 2(-B ,b), 
which gives the claimed bound using (14.51) . 

The second term is estimated using directly Cauchy-Schwarz estimate and the fact that 
sup [0jb] \L'(x)\ < CW~\ We get 



\B k (v)\ <C-W 1 ||G r fc||i,a(o J 6) • ||t/||l,a(-B ,6). 



That gives the claimed bound using again (j4.5p . 



□ 
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Define F := lyyd x u and G := ly/d y u. By definition, F k {x) is the Fourier coefficient of 
the function F(x, ■) with respect to the Fourier basis t— > cos yk^j^J^ 

Using Plancherel formula we get 

^F fc (x) 2 ^< [ L(X) \F(x,y)\ 2 dy. 



k>l 

For the same reason, but using this time the sin basis we have 

,L(x) f L{x) 



T,(r\ f^ x> 

k>i ^° 



Integrating with respect to x and bounding y from above and L(x) from below uniformly 
we get the following lemma : 

Lemma 4.2. For any bo there exists C depending only on the billiard and bo such that, for 
any b < b , 

(4-9) E ll^lliw ^ll^lliw 

fc>i 



(4-10) 

1,2(0,6) — C\\d y ii 



2 

y a \\L2{w b )- 

fc>i 



We now switch to the control estimate. We begin by dealing with the modes for which 
^f-E>E. 

4.3. Large modes. 

4.3.1. A control estimate. Equation (14.31) may be rewritten 

(4.H) - Ufc + [j2^T ~ E J U k = h k , 

where hk is the element of H~ l defined by 

(4.12) h k :=A h + ju' k 

The H^ 1 norm of h k is now estimated using the following 
Lemma 4.3. There exists a constant C := C(b ) such that for any b < bo and any k such 

1,2 2 

that ^-E>E the followinq estimate holds: 

(4.13) ||Mff-i(_B ,b) - C ( b o) ( kvy \\ F k\\v(o,b) + &7_1 11^11^(0,6) + & 7-1 |KIU2(o,&)) . 
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Proof. Using remark 14.21 the norm of Ak is uniformly controlled by the norm of and 
the latter is estimated using Lemma (14. ip . To estimate the if -1 norm of j;u' k , we first set 



v = and remark that 



v 



V 



T 



(L 



l\2 



L 2 



We choose a test function and estimate 

,6 /£/ 



h 



B 



L 



Uk I 4>dx. 



We perform an integration by parts, use that ^^j- < Cb" 1 1 l x >o, then apply the Cauchy- 
Schwarz inequality to get 



Bo 



Uk) 4>dx 



L'(x) 



-B L ( X ) 



UklX) 



b 1 (x) dx 



We then estimate 



Wfc || L 2 (O,&)||0 \\L 2 (-B ,b)- 

h '!?{?) (L'(x)) 2s 



Bo 



L(x) L 2 {x) 



u(x)(j)(x) dx 



We perform an integration by parts, use that 

'L"(x) (L'(x)) 2 
L(x) L 2 {x) 

then twice apply the Cauchy-Schwarz inequality to get 



< Cxl~ 2 , 



h C 



-Bo \J0 

< Cb 7 ' 1 



a~ 2 Koi^) 

II L 2 (0,6) 1 10' || L2 (0,6)- 



The claim follows using (14. 5p . 



□ 



The variational formulation of equation ( 14. lip is given by 



(4.14) 



u' k v' dx + 



Bo 



2^2 



k 2 7T 

L 2 (x) 



— E UkV dx = hk(v). 



i 2 2 

Since ^M? E > E, the left-hand side of the preceding equation is a continuous qua- 

L o 

dratic form on H^—Bo^b) so that, using Lax-Milgram theory, there exists a unique Vk in 
Hq(— B ,b) that satisfies (14. lip in the distributional sense. 

The following lemma allows us to estimate the L? norm of this Vk- 
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Lemma 4.4. There exists a constant C depending only on bo but neither on b < bo, k, nor 

E such that, if E > 1 and ^ E > E then the variational solution v k in Ho(—Bo,b) to 

equation (4-H ) satisfies 

(4.15) KIU 2( -Bo,&) < C(bo)(^\\F k \\ LH o^+E~h^G4 LH o,b) 

+ J ET^ 7 ~ 1 |K|| £2(0,6) 



Proof. Since v k is a variational solution, putting v = v k in (14. 14ft we get 
(4.16) [ \v' k (x)\ 2 dx + [ (^^-E^\v k (x)\ 2 dx = h k (v k ). 



B J -B 



In the regime we are considering the second integral on the left is positive so that we obtain 

\v k (x)\ 2 dx < \h k (v k )\ < \\h k \\ H -i\\v k \\ H i. 



-Bo 

Since v k is in Hq(—Bo, b), Poincare inequality gives c(b) a positive continuous function of b 
such that 

b 

\v k (x)\ 2 dx > c{b)\\v k \\ 2 m . 

-B 

This gives a constant C depending only on 6 such that, for any < b < b , we have 

IKIU 1 < C\\h k \\ H -i. 
We now use again equation (14.161) to obtain 

\ rb 

y-j2 E J J \v k (x)\ 2 dx < \\h k \\ H -i\\v k \\ H i < C\\h k \\ 2 H -i 

with the preceding bound. Using estimate ( 14.13P we obtain 

(-j2 E J lkfelU 2 (-B ,6) < C (kW\\F k \\ L 2 {0ib ) + b^ 1 \\G k \\ L 2 {0 ,b) +b 1 ~ 1 \\u k \\ L 2 m ) . 

f 1,2 2 \ 2 

We divide both sides by ( ^ E 1 . The coefficient in front of b 1 ||i r fc||£2(o,&) is bounded 

by using the fact that 

k 2 L 2 ( E 

sup w = sup — 1 + — 

k ^~E>E ^l\~ E Z > E 71 V ^ 



7T 



2 



1 , 

E 



For the two other terms, we use simply that - E > E. This gives the lemma. □ 
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We can now let Wk = — v^. By construction, Wk is a solution to the homogeneous 
equation 

(4.17) _W + (*g- B )„ = 0. 

Moreover, since both and Vk satisfy Dirichlet boundary condition at — B we have 
that Wk(—B ) = 0. 

Since the 'potential' part in equation (14.1 7p is bounded below by E, concentration prop- 
erties of solutions may be obtained using convexity estimates. 

Lemma 4.5. For any b < bo, any solution w to h4-l r ty such that w(—B ) = satisfies 

/b pb 
\w\ 2 (x) dx < (B + bo) / \w\ 2 (x)dx. 
■B JO 

Proof. Multiplying the equation by w we find 

' fcV 
l?{x 

It follows that 

{w 2 )" > f3 2 w 2 , 

for some positive /3 (here /3 2 = 2E). 

Since w(—Bq) = 0, using the maximum principle on [— Bq,£\, we obtain for all —Bo < 

x < £ < b 

2 2 sinh(/3(x + B )) 



k 2 n 2 

w"w + ( — r-- - £ ) w 2 = 0. 



w J (x) < u/(£)- 



sinh(/3(e + i?o)) 

For any t G [0, 1], define x(t) = —Bq + £(-Bo + b) and = tb. Since for any £ we have 
—Bo < x(t) < £(£) < &0j w e may integrate the preceding relation : 

9/ / , f 1 9/*/ sinh(/3(x(t) + 5 )) , 



o Jo 



Since sinh is increasing the quotient of sinh is bounded above by 1 and we obtain 

/b pb 
w 2 (x) dx < (Bo + b) w 2 (x) dx. 
-B JO 



□ 



Putting these two lemmas together we obtain the following 
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Proposition 4.6. There exists a constant C depending only on bo such that for any b < bo, 
for any k and E such that k 2 ir 2 jL\ — E > E and E > 1 

(4.18) \\u k \\ L 2 { _ Boib) < C (b^\\F k \\ L 2 m + E~h^\\G k \\ L 2 m 

+b~^\\uk\\ L 2(o,b)) 

forC = C(b ). 

Proof. According to Lemma 14.51 we have 

||wfclU 2 (-Bo,&) ^ Cb~2\\w k \\ L 2( Qjb) , 

where w k = u k — v k and v k is the variational solution constructed above. Using the reverse 
triangle inequality, we obtain 

_ i i _ i 

\\u k \\ L 2(-B ,b) < Cb 2 |KHl2(o,&) + (C + 6 2 )6 2 \\v k \\ L 2(- Bo ,b)- 

The claim will follow using estimate f 14 . 1 5 1) of Lemma 14.41 Observe that the prefactor of 
||wfc||L 2 (-s ,6) * s a ^ fi rs t ( U P to a constant prefactor) 

6-s + b-^E-h^- 1 . 

Since E'^b 1 " 1 is uniformly bounded we obtain the given estimate. □ 

4.3.2. Summing over k . We will now sum the preceding estimates over k. We thus introduce 

( kwy\ 
\L(x)J 



u + (x,y) = ^2 u k (x)sm 



*Z2!—E>E 
l<\ — 



and prove the following proposition. 

Proposition 4.7. There exists bo and E and a constant C depending only on E and b 
such that if u is an eigenfunction with energy E > E$ and b < bo, then : 

\Mh {R) < c [(ft^HMIwi + E- l b 2 ^\\d y u\\l 2[Wb) + b-'WuWl^ 

Proof. We square estimate f l4.18p . sum over k, and use (13. 2p and Lemma [4.21 □ 

Observe that the controlling term in the preceding estimate is supported in the wing 
region. However, compared to the usual bounds (such as in [6]) there is a loss of derivatives 
since we need d x u and d y u in the wings. 

We also obtain the following corollary. 

Corollary 4.8. Let bo and Eq be fixed, there exists C depending on the billiard bo and Eq 
but not on the eigenfunction nor on b < bo such that 



\ u +\\l 2 (r) ^ C 



(b 2 ^E + b^)\\u\\l Hn) +6- 1 N|£ a(H0 
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Proof. We bound ll^^ll^^) an d H^^lli 2 ^) ^ -^ll M lli 2 (^) an< ^ use ^ ne ^ ac ^ ^ na t the norm 
over Wb is less than the norm over W. □ 

It remains to choose b in a clever way to obtain the desired bound. 

4.3.3. Optimizing b. We will choose b to be of the form M~ 1 E~ a for some constants M 
and a to be choosen. As long as a is positive, there is some large Eq such that for any 
E > E then b = ME~ a < b so that we can use the preceding proposition. 

We obtain 

(4-19) \\u + \\h (R) < C [(M 1 - 2 ^ 1 -^- 1 ) + E-^-V)\\u\\l HQ) + ME«\\u\\h {w) 



It remains to make good choices to obtain the following proposition. 

Proposition 4.9. There exists E and C depending only on the billiard such that for any 
u eigenfunction with energy E > Eq the following holds : 

(4-20) IMl!a(fl) < \\H\ 2 L 2 [n) + ce^\\u\\I hw) 

Proof. We choose a := and M such that CM 1-27 = |. For £ large enough, E~ a ( 2 ~'~^ 
goes to zero. It is thus bounded by for E large enough. Replacing in (I4.19P we get : 

\\u + \\ 2 L 2 {R) < ^||u||^ (n) + CE^\\u\\ 2 L2{w 



The claim follows. □ 



4.4. Small modes. We now consider modes for which ^ E < E, and this time we 

rewrite the equation P k {u k ) = A k in the following form: 
(4.21) - u k - z k u k = h k , 

,2 2 

in which we have set z k := E ^- and 

a L' , k 2 ( 1 1 \ 
h k := A k + —u k + — ( - — j u k . 



7T 



4.4.1. The control estimate. Since z k > —E we can use the results of the appendix to 
control |K||L2(-B ,f>). 

To do so, we need to estimate the norm of hk in H^ 1 {—Bq 1 b). 

Lemma 4.10. There exists some constant C depending only on b such that, for any b < b 

1,2 2 

and any k such that E < E, the following holds : 

(4.22) ||Mh-i(-b ,6) ^ C(kW\\F k \\ L 2 m + ^- 1 ||G fc || i2(0ib) 

+ (V- 1 + k 2 W +l )\\u k \\ LHm ) . 
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Proof. We use the definition of hk = j^-^k + j;u' k + ^ — -pj ut and estimate each 

term separately. The first term is estimated using (14. ip and remark 14.21 The second term 
is estimated as in the proof of lemma 14.31 

The same method applies to estimate the third term. We introduce 



Uk(x)4>(x) dx 



B \Lq L 2 (x) 



and observe that {j?— pr^y ) is 0(x\). Integrating by parts and using twice Cauchy-Schwarz 
inequality gives 

J 3 < C6 7+1 ||0 / || jL 2( O) b)||«fc||i2(o ) 6). 

Using the definition of the if" 1 norm (see (14. 5ft ) and putting these estimates together 
yield the lemma. □ 



For any Bel, denote by v(E) := min j 
Remark 4.3. Observe that v(E) < min j 



E - 
E 



k 2 n 2 


l 2 n 2 


T 2 




fc 2 7T 2 


7T 2 







have 
(4.23) 

for some constant c. 



, (M)eNxN 

, k G N j , so for E large, we 



v(E) < cVE 



Lemma 4.11. For any (3 > 0, there exists some c such that the following holds. For any 



k such that z k = E — -M- > (5 



u(E) 



sin(E ox /i^)| > c 



Proof. First we use that there exists some c such that 

Va; G K., |sinx | > c • dist(x, 7rZ). 
We denote by Ik the integer such that 

dist f y/zk~, = v^fe ~ 
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so that we have 



'Zk - 



ZfcTT 

So 



> C- 



Zk ^ Bo 



E 



> c- 



k 2 n 2 _ It* 1 



where, for the last bound, we have used the lemma 14.121 below. 
The claim follows by definition of v{E). 



□ 



Lemma 4.12. Fix a > and denote by I the (step-like) function on [0, oo) that is defined 
by 

|A — l(X)a\ = dist(A, aZ) 
Then there exists some C such that 

VA G [0,oo), X + l(X)a < CX. 

Proof. Define / by /(A) = A+ ^ A ^" . First, since I vanishes on [0, |], we have /(A) = 1 on 
this interval. Second, the function / tends to the limit 2 when A goes to infinity. Finally, 
on [f , M] we have /(A) = 1 + l -^a < 1 + ^Jtl. □ 



Putting these estimates together, we get the following 

Proposition 4.13. There exists bo and Eq and a constant C := C(bo,Eo) such that the 

following holds. For any E > Eq, for any k such that E < E and for any b < b , we 

have the following estimate 



(4.24) 



\Uk\\L 2 (-B ,b) 



< C 



e\ 



v{E) 



EW + 2\\F k \\ L 2 m +V-i\\G 



r fe||i 2 (0,6) 



(1 + EW +2 )b-^\\u k \\ L i 



(0,6) 



Proof. For any k we let z k = E 



k 2 W 2 



and use the estimates of the appendix combined 



with the bound on h k given by Lemma I4.1UI For k such that z k corresponds to estimates 
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\Uk\\L 2 {-B ,b) < C 
< C 



62 \\h 



k\\H~ 1 {-B ,b) 



+ b 2 ||'U fc || L 2( 0ifc ) 



/c6 7+ 2 ||F fe || L 2( ,6) + b 1 2 ||G fc || i2(0ib) 



+ (ft7 + ^7+2 + ^ H|K|| L 2 



(0,6) 



We now use that k = 0{E?) in the regime we are considering. We also remark that 
{W + k 2 W +2 + 1) = 0(1 + EW +2 ). 

Otherwise (i.e. for k such that zj, corresponds to estimate (lA.llj) ). we have to add a 
global I sin(5oA/^fc)| _1 prefactor. Using Lemma 14.1 1[ we have 



sin^o^r 1 < C^- < C E 



u{E) ~ is(E) 
We thus obtain that for any k the following holds : 

E\ 



KI|z,2 ( _,B 0ife) < C ■ max 1, 



u{E) 



En^\\F k \\ L2 



(0,6) 



+b^\\G k \W m + (1 + EW +2 ) b-*\\u\\ L 2 m 
Using (I4.23p . for large E we have is bounded from below, so that the claim follows. 



□ 



4.5. Summing over k. We use the estimates of the preceding sections to obtain a control 



on llit-llra/m in which we have set 



ilnl-E<E 



We prove the following proposition. 



Proposition 4.14. There exists bo and E and a constant C depending only on E and bo 
such that if u is an eigenf unction with energy E > E and b < bo, then 



\U- 



\L*(R) 



< C 



E 



^ +1 ||^||| 2(w) + 6 2 ^ 1 ||^||| 2(w , ) 



v(E) 2 

+ {i + ew+ 2 ) 2 b^\\u\\ 2 LHW) 
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Proof. We square (14.241) and sum with respect to k. The Lemma l4~2l controls ^2 ||-^fc|| 2 and 
^2 ||Gfc|| 2 - Plancherel formula takes care of ll M A;|| 2 - We also use as before that the norm 
over Wb is smaller than the norm over W. □ 

As for the large mode case, we get a corollary using the fact that ||9 x m|| 2 and ||9j/u|| 2 are 
bounded above by i?||tt|| 2 2 ^. 

Corollary 4.15. There exists b Q and E Q and a constant C depending only on E and b 
such that if u is an eigenfunction with energy E > Eq and b < bo, then 



\u. 



\L*(R) 



< C 



E z 



HE) 



E 2 



HE) 



u\ 



L 2 (Q) 



(l + ^ +2 ) : 



Eb- 



HE) 



2 \\ U \\l^{W) 



4.6. A non-resonance condition. We now want to make the previous estimates explicit 
with respect to E and b so that we can use a similar optimization procedure as for the large 
modes case. We thus impose some condition on v(E). Namely, for any e > 0, we introduce 
the following set 



Z £ := {E e R,\v(E) > c E- £ } 



E e R | Vfc, I e N, \E- — 2 -2-| > c E~ 

In other words, the set Z £ consists in energies that are far from the Dirichlet spectrum 
of the rectangle [—-Bo, 0] x [0, L Q }. It is natural to say that such energies are not resonating 
with the rectangle. The coefficient Co which is irrelevent when e > has been chosen in 
such a way that Weyl's law for the rectangle implies that Z is not empty. Note however 
that, although expected, it is not clear that there actually are eigenvalues in Zq, nor for 
that matter in Z £ . 

Once e is fixed, the estimate of the corollary 14.151 becomes : 

h-Whm < c ■ [ (b^E^ + b^E^) \\ u \\l Ha) 

(4 ' 25) +(l + Eb^) 2 b-'E^\\u\\l Hw) 

4.6.1. Optimizing b. As before we let b = ME~ a for some positive a and try to optimize 
the bound. 

Proposition 4.16. Define a by 

'3 + 2s 2 + 2e 



a = max 



2 7 + l ' 2 7 - 1 
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There exists E Q and C such that for any u eigenfunction with energy E in Z e such that 
E > Eq, the following holds : 

(4.26) \\u-\\l* {R) < \h\\h<n + c-e 1+ * + °- 

Proof. With the given choice of a it is possible to choose M so that the prefactor of ||u||f^(™ 
is i for E large enough. The claim follows remarking that the definition of a implies 
a > j^TT > ~To so ^ na ^ the prefactor (1 + Eb 1+2 ) 2 is uniformly bounded above. □ 



5. Non-concentration Estimate 

We now put all the estimates together to obtain the following theorem. 

Theorem 2. Fix e, and define p by 

/2 + 7 + 2( 7 + l)e l + 2 7 + 4 7 e 

p := max , 

y V 2 7 + l '47-2 

There exists Eq and C such that any eigenfunction u of Q with energy E in Z E such that 
E > E satisfies : 

NIl^o) < C ■ E p \\u\\ L 2 {w) . 

Proof. We first remark that whatever the exponent a is we always have l+2e + a>l> 
2^ry so that the exponent for the small modes is always larger than the exponent for the 
large modes. Thus, adding the estimates from propositions 14.91 and 14.161 we obtain 

htvm ^ \\\ u Wl^) + CEl+2e+a \Hh(wy 

Since \\uf L2(m = ||u||2 2 - we get 



\u\\ 2 ma) <(i + CE 1+2e+a )\\u\* 2 



\l 2 {r) — W u \\l 2 (U) H^llL^vy) 
1 

2 ll "Hi 2 (0) ^ \ L 1" )\\ U \\L 2 (W) 

When E is large the constant 1 can be absorbed in the term with a power of E. The claim 
follows by computing 1 + 2e + a for both possible choices of a and taking square roots. □ 



We state as a corollary the corresponding statement for the Bunimovich billiard (see 
theorem [T]) . 

Corollary 5.1. In the Bunimovich stadium, for any e > there exists Eq and C such that 
if u is an eigenfunction of energy E in Z E such that E > Eq then the following estimate 
holds: 

IM|z,2 (n) < CE e || M || i2(iy) . 

Proof. We let 7 = 2 so that a = max(^ £ , ^p). For any non-negative e, -^p < this 
makes the proof complete. □ 
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Remark 5.1. The bounds in [6] gives a similar control with 1 as the exponent of E. 
Our bound thus gives a better estimate as long as s < |. As it has been recalled in 
the introduction, it is quite natural that the non-resonance condition allows to get better 
bounds. 

Remark 5.2. We could deal with quasimodes by adding an error term to A that is con- 
trolled by some negative power of E. There will be mainly two differences in the analysis. 
First the second term A will not have support away from the rectangle anymore and second, 
in the optimization process, we will have to take care of the new error term (which will 
possibly change the range of applicable exponents). 

Remark 5.3. By adding the estimates in propositions 14.71 and 14.141 we get a different 
control estimate, where the control still is in the wings but now with a loss in derivatives. 
We haven't tried to optimize this bound. 

Appendix A. One dimensional Control Estimates 

The aim of this appendix is to provide a control estimate for the equation 

—u" — z ■ u = h 

on [— B , b] of the following type 

||'"||i2(_ J Bo,0) < Ci\\h\\ H -i(-B ,b) + C2IMI 1,2(0,6), 

in which we want an explicit dependence of the constants G\ and C*2 on z and b. It is now 
standard (see [7]) that if b is fixed then we can choose C\ and C2 to be independent of z 
but what we need is an estimate when b goes to 0. 

We first need a few preparatory lemmas. 

Lemma A.l. For any e > 0, there exists a constant C := C(e) such that for any b, for 

any h £ H^ 1 {—Bq, b) and any z such that z < there exists a solution v p £ -£^)(0, b) 

to 

-v'p - zv p = h, 

in T>'(0, b) and 

(A- 1 ) \\v P \\mo,b) < Cb\\h\\ H -n-Bo,b)- 

Proof. First we note that h, when restricted to (0,6) also belongs to i7 _1 (0,6) and that 
II ^11 h- 1 (o,b) ^ II^H-ff-if-Bcb)- The proof follows from a standard resolvent estimate since, 
on (0,6), the bottom of the spectrum of the self-adjoint operator v 1— > —v" with Dirichlet 

2 

boundary condition is p-. We include it for the convenience of the reader. We decompose 
v p in Fourier series : 

v p( x ) = ^Qfc sin(-^x). 

k>l 



We have 



hence 
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k 2 7T 2 \ ,kll 



fc>i v 



b 2 y b 



2 _ V V b " J | 2| 

H- l (0,6) - fe 2 7r 2 



fe>l I 3 " 



or 



2 \ ^2i-2 

H-i(0,6) ^ n 



'l-#2) 2 
fc 2 7T 2 / 



inf 

k>i ir 2 k 2 



\v II 2 



\ 2l-2|| ||2 

> cix b \\v p \\ L 2 { ^ b) . 

The claim follows since the inf is bounded away from zero in the regime we are considering. 

□ 

Lemma A. 2. For any z < ^~^ n Let w G Hq(—B ,o) be a solution to 

II n 

— w — zw = (J 

in T>' {{— Bo,b) \ {0}) . Then there exists a constant A such that w = AG, in which the 
function G is defined by: 

sin( v / i(a:+Bo)) sin( v / 2fe) 



(A.2) G(x) 



—r 7= — , x < 0, 

V2 



sin(y/z(b-x)) sin( N /zBo) g 
y/z s/z ' 



Proof. Let u> be such a function then necessarily there exist two constants A± such that 



sin( v / i(a:+Bo)) x < Q 



By assumption u; G -ff 1 and hence is continuous at 0, so 

A sin(y^Bo) _ ^ sin(v/i6) 



'2 

In the regime we are considering sm (v^ b ) ^ o, hence we can divide by this expression and 
express A_ in terms of A + . The claim follows. □ 

We finish these preparatory lemmas by establishing the control estimate for multiples of 

G. 

Lemma A. 3. 
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(1) For (3 such that < (3 < there exists B\ = Bi(f3) and C := C(/3) such that for 
any z < j3 2 and any b < B\ the following estimate holds 

(A.3) WGWvi-Boja) < Cb-*\\G\\ 

(2) For any /3,e > there exists B\ := B\{0) and C := C(/3,e) such that, for any 

b 2 



b < B\ and (3 2 < z < - the following estimate holds : 



b' 2 

( A - 4 ) I|G'||l2(-Bo,0) < C , ||<2||z,2( 0& ). 

sin(VzBo) 



Proof. In case (1), we first assume that z < — Zq for some positive Z . we set z = —u 2 and 
compute 



/° \G{x)\ 2 dx = T sinh 2 (a;(So + x))dx 

J-Bn U J-Ba 



Bo 







|Gr(a;)| a da; = \ — — / smh. 2 (u(b — x)) dx, 

0J Z 



By a straightforward change of variables we get 

-0 n;„l,2/ ( .,L\ /-wSo 



/ G x dx = ^ / smh 2 (Ode, 

J-B W ^0 



Bo 

2 



sinhV^o) 
F 



|G(x)Mda; = / sinh 2 (0df. 

o ^ Jo 



We set := ^ J^x) — so that finally, we obtain : 

\G{x)\ 2 dx = F( ^ B ^ ■ f \G{x)\ 2 dx. 
B F(ub) J 

It is straightforward that F(X) is positive, tends to 1 at infinity and that F(X)/X tends 
to | at 0. As a consequence, there exists some C(Zq) such that, for any z < —Z 2 , 

\G{x)\ 2 dx < Cmax(l,(u;6)- 1 ) / \G(x)\ 2 dx, 



For b < B\ and u> > Z , we have max(l, (cob) < max(l,u; v )b 1 < Cb 1 which gives the 
claim for this range of parameters. 
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We now assume that we have — Z% < z < f3 2 . We have 



\G(x)\ 2 dx 



B 



sin(y / i6) 


2 M 


r 

sin(y^(a; + B )) 




J-B 





= b 2 

Cb 2 



sm(y/zb) 



J~zb 



2 ^0 



So 



sin( x /i(x + B )) 



dx 

2 



v/i(x + B ) 



(x + B ) 2 dx 



where the constant C comes from the fact that the function sin ( w ) i§ continuous and its 

w 

argument belongs to a fixed compact set. On the other hand, by a simple change or 
variables we have 



\G(x)\ 2 dx 



sin(- v /i5 ) 


2 »b 


sm(y/zx) 




Jo 





dx 



,b 3 

> cB 2 -, 



in which c is given by 



sm(y/zB 



inf 

0<x<B! 



sin(\/zx) 



y/zX 



Using that sm ^ is continuous and does not vanish on (— oo, 7r) and choosing Bi accordingly 
we obtain the first bound. 



For case (2), we first use homogeneity and prove the bound for G := zG. We have 



So 



G(x) dx = |sin(v / i6)| / |sin(\/i(a; + -B ))| dx 

J-B 
2 



< B \sm(X)\ 

in which we have set X := \fzb. On the other hand we have 

G(x) dx = \sm(y/zB )\ / |sin(\/J^)| dx 

Jo 

x 



6-|sin(v^i? )| 2 -l^ |sin(0| 2 ^ 



with the same X. Under the assumptions, X belongs to a compact subinterval of [0, 7r). 

x\ S in(x)\ 2 Jo* I s ^ n (0l 2 ^ * s continuous, the claim 

□ 



Since on this interval the function X y 
follows. 



We can now prove the following 
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Proposition A. 4. There exist (3 and Bi := Bi(/3), such that ifb < B\ andv G Hq(—Bo,V) 
satisfies 

—v" — zv = h, 

with h that vanishes on (—B Q ,0), then the following estimates hold: 
(1) If z< (5 2 , then 



(A.5) \\v\\ L H-B fi) < Ci 

(2) If(3 2 <z<±, 

(A.6) \\v\\m-B ,o) < Ci 

(3) If^<z then 



6i 



sin(5 0A /i)| 



H-i(~B (h b) + 6 2 IN 



|^||h-1(-B 0) 6) + 



L 2 (0,6) 



6-1 



sin(5 -v/i)| 



lklU 2 (o,b) 



(A.7) 



I^||l 2 (-b ,o) 



< C, 



& 2 llfc||ff-i(-B ,6) + & 2 IMIl 2 



L 2 (0,b) 



Proof. In the first two cases, we have z < h < ?j. We may thus consider v p as given by 
lemma |A~T1 and define v p by extending v p by for negative x. Observe that w := v — v p is 
in Hl(—B 0) b) and satisfies 

— io" — zu; = 

in T>((—Bo,b) \ {0}) so that u — 5 P = AG for some A according to Lemma IA.21 Using 
Lemma IA.3I we obtain in the first case 

Ik - V p \\tf(-Bo,0) < Cb~~*\\v - Up ||ia(o } 6). 

We use the triangle inequality on the right-hand side and the fact that v p is for negative 
x and coincide with v p for positive v. We obtain 

\\v \\m-b ,o) < Cb~^ (|H| L 2( 0)6 ) + \\v p \\ L 2 m ) . 

The claim then follows from the estimate on || ^ p || i 2 (o,6) i n lemma |A~T1 We prove the second 
case by following the same argument, inserting the corresponding bound for G. 

The third case will follow the same lines but we will introduce a different particular 
solution v p , following then even more closely the proof of [TJ. We set A = yfz. 

Denote by H the unique L 2 function on (—Bo, b) that vanishes on (—Bo, 0) and such that 
H' = h in the distributional sense. The L 2 norm of H is related to the B~ x norm of h by 
the relation 



\H 



H(y)dy 



\L 2 (-B ,b) 



\h\ 



H-\-B ,b)- 



The Cauchy-Schwarz inequality then implies that 



(A.8) 



\H\ 



L 2 (-B ,b) 



> (l + b 



\H-l(-B ,b)- 
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Set 



•/-Bo A 



then v p satisfies 



// \ 2 ttI 

— 1> — X v„ — H 



in V'(—B , b) and v p (—B ) = but the boundary condition need not be satisfied at b. We 
thus have 

v(x) = v(x)-v (b f™(\(x + B )) 
v{x) v p[ x) v pW sin ^ Bo + b) y 

The function v — v p is thus a multiple of sin(A(x + -Bo))- 
We have 

r-0 



|sin(A(x + J Bo))| 2 dx < B 



Bo 



and 

J \sm(\(x + B ))\ 2 dx>±(b-±). 

Hence, in the regime under consideration we have 

(A.9) \\v - u p ||i3(-fl 0l 6) < Cb~^\\v - v p \\ L 2 { Q fi) . 

We perform an integration by parts in v p and observe that the boundary contributions 
vanish because H vanishes near —Bq and sin(A(|/ — x)) vanishes at y — x. 

Finally, we obtain 



v p (x) = / cos(X(x - y))H{y) dy. 

J-Bo 

It follows that v p is identically on (—Bq,0) and that, on (0,6), it satisfies 

\v P {x)\ < \\H\\ L 2 { _ Bo>b) ^x~. 
Squaring and integrating, we get 

|kp||L 2 (0,b) < b\\H\\ L 2(_ Bo ^. 

Using the triangle inequality in (1A.9|) and inserting this bound, the result follows for b < | 
using (IA.8p . □ 

In the paper, we will need to relax the condition that v(b) = 0. This can be done using 
a standard construction related to a commutator method. We will get the following 
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Theorem 3. There exist (3 and four constants B\, C±, C2, C3 depending only on (3 such 
that the following holds. For any b < B\, for any function u in H x {— Bq, b) that satisfies 

—u" — zu = h, 

with h G H~ l (—B ,b) and such that u(—B ) = and h vanishes on (—B o ,0). Then, the 
following estimates hold: 

(1) Ifz < f3 2 , then 

(A.10) IMUvbo.o) < Ci 

(2) Iff3 2 <z<±, then 



b 2 \\h\\ H -i(-B ,b) + b M u IU 2 (o,b) 



b\ 



b~\ 



\MB^\ mH -^ + \MB ^)\ ML2m 



(A.ll) \\u\\ L 2 { -b ,0) < C x 

(3) If h < z, then 
(A.12) ||u||i2(-Bo,o) < C 3 

Proof. Define a smooth cutoff function pi such that pi(x) is identically 1 if x < | and 

- p;,u then v G 



6*1 



|h-1(-So,6) + 6 2 IMU 2 



(0,6) 



J I UUJ . J-JClllie a SlilUULll UULUll 1U11CL1U11 SUU1 

identically if x > 1 and let pb be the function x 1— >■ Pi(f)- Define t 
//q(— i?o,&) and satisfies 



—v 



■ zv = h + 2(p' b u)' - p" h u. 

The right-hand side vanishes on (—-Bo, 0) so that, in order to use proposition IA.44 we have 
to estimate its H" 1 norm. The strategy is the same as in the proofs of lemmas 14.31 and 

EM 



LUJ 

An integration by parts followed by the use of the Cauchy-Schwarz inequality 



gives 



-B 



< 



\p'u\\L^,b)\\<t>'\W 



b 



Thus, 



C, 



IKpWIU-i < T-|hlU2( ,6)- 



' 

The third term can be estimated using the same method. Indeed 



p'b u 4> 



'0 \Jo 



< 



p'b(y) u (y)dyj <P'( x ) dx 
pb(y) u (y) d y\\L2(o,b)H'\\L2- 
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Using again Cauchy-Schwarz inequality and the fact that \p'l(y)\ < Cb 2 we get 



p'b(y) u (y) d y 



< Cb 2 \\u\\ L 2( ^y/x 



We obtain 



pb(y) u (y) d y\\^(o,b) < cb 2 \\u\\ L 2 i0jb) \\^/x\\ L 2 m 

< Cb' l \\u\\ L 2 m . 



It follows that 



\\h + 2(p' b u)' - p'bU\\ H -i(-Bo,b) < ||^||h-i(-Bo,6) + Cb 1 ||«||£2(o,&)- 

We obtain the theorem by plugging this bound into the estimates of the proposition [AT4l D 
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